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I. INTRODUCTION 

A characteristic property of the majority of quantum field theories (QFT) is the divergence 
of integrals appearing in the perturbation theory for the calculation of physical quantities such 
as mass and charge of the interacting particles. The divergences appear in the integrations over 
momenta in intermediate states both on the lower limit (infrared divergence) and on the upper 
limit (ultraviolet divergence). In order to circumvent this difficulty, a renormalization proce¬ 
dure is used, which allows the redefinition of the initial parameters of the system through their 
observable values. The renormalization scheme was firstly developed for quantum electrody¬ 
namics (QED) [1] and later generalized to other QFT models [2], These schemes can be used for 
so-called renormalizable theories, for which the reconstructed perturbation theory can be built 
in a way that the infinite values are included in the definition of “physical” charge and mass 
and, therefore, do not appear in other observables of the system. However, even the founders 
of QED anticipated “that the renormalization theory is simply a way to sweep the difficulties 
of the divergences in electrodynamics under the mg.” (R. R Feynman [3]). In many papers P. 
A. M. Dirac wrote that this approach was in contradiction with logical principles of quantum 
mechanics [4, 5], 

Accordingly, the question arises whether these divergences are an intrinsic property of quan¬ 
tum field models or they are caused by the application of perturbation theory for the calculation 
of physical quantities, which are non-analytical functions of the coupling constant such as in the 
theory of superconductivity [6]. A large number of works is devoted to this problem, neverthe¬ 
less a solution still has not been found up to now. However, this question is of great importance 
for a correct mathematical formulation of fundamental physical theories and is essential for ex¬ 
amining the applicability of non-renormalizable theories [7] for the description of real physical 
systems. 

Let us recall that in standard perturbation theory the Hamiltonian of non-interacting fields 
is used as a zeroth-order approximation, while Fock states of the free fields are employed for 
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the calculation of the transition matrix elements in the subsequent corrections to observable 
characteristics of the system. This approach is based on the assumption of an asymptotic switch 
off of the interaction between the fields [8]. However, in a series of works it has been shown [9- 
12] that the infrared divergence arises just because of the use of asymptotically free field states. 
As follows from reference [9] and the subsequent publication [10], the infrared divergence 
disappears in all orders of perturbation theory in QED if, in the zeroth-order approximation, the 
coherent states of the electromagnetic field bound to the particle are used and the parameters of 
these states are appropriately chosen. 

At first glance, this may contradict representation theory in quantum mechanics, in accor¬ 
dance to which the result of the calculation of the observable characteristics of the system should 
not depend on the choice of basis states, provided that those form a full basis in a Hilbert space as 
for free field states. However, this statement is correct only for the exact solution of the problem, 
whereas individual terms of the perturbation series can change with a different choice of basis 
in zeroth-order approximation. As was demonstrated in reference [13, 14] the transition from 
one basis to another corresponds to the partial summation of a divergent series within standard 
perturbation theory and allows for the non-perturbative calculation of subsequent corrections 
in the form of a convergent sequence. A good example of how the basis choice influences the 
approximate calculation of the characteristics of the quantum system with a continuous spec¬ 
trum is given by the scattering at a Coulomb potential [15]. In this well known case, the wave 
function of the system has no singularities. In contrast, Born’s scattering amplitude, approxi¬ 
mately calculated via an asymptotically free basis, displays a singularity for scattering at small 
angles. As was demonstrated in references [15], this singularity does not appear with the use of 
non-asymptotic wave functions. 

The main goal of our work is to investigate whether a proper choice of basis in zeroth-order 
approximation allows to constmct a calculation scheme free of ultraviolet divergences. In order 
not to overload the proposed approach with details related to the internal degrees of freedom and 
to render all calculations as transparent as possible, we investigate as a representative example a 
model system, which consists of a non-relativistic particle without spin interacting with a scalar 
quantum field. A standard-perturbation-theory series, in this case, does not exhibit infrared 
divergence, contains however ultraviolet divergence. This results in a dependence of the energy 
of the ground state on the undefined momentum cut-off, which is required for the calculation of 
high-order corrections. Consequently, our task is to prove that the energy of the ground state 
of the considered model system can be calculated without any additional parameters such as a 
momentum cut-off. At the same time, it is important to show that the energy of the system is a 
non-analytical function of the coupling constant and consequently can not be represented as a 
series in the framework of conventional perturbation theory. 

With the inclusion of a field polarization our employed model coincides with non-relativistic 
QED [16] or if the field is scalar it has a physical realization in solids [17], where however, due 
to the discrete structure of a crystal, a natural momentum cut-off intrinsically appears, defined 
via the Brillouin-zone boundary. In free space this regularization is not present and has to be 
artificially included, e.g., via lattice models [18], where the boundary momentum is defined 
through an artificial lattice period. In contrast, in our formulation we will consider a system in 
free space without neither natural nor artificial cut-off. 

In addition, in a series of works [19-22] an analogous model of a particle interacting with 
a scalar quantum field with the momentum cut-off was used for the investigation of the funda¬ 
mental mathematical problem of the existence of the solutions of the Schrodinger equation. 

The article is organized in the following way. In section II the model of a non-relativistic 
particle without spin interacting with a scalar quantum field is described and its parameters 
are calculated in the framework of conventional perturbation theory. In section III the basis of 
non-asymptotic states is investigated and the iteration scheme of the calculations is presented. 
The zeroth-order approximation, which is found to be free of ultraviolet divergence for the 
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energy and effective mass is then worked out using this basis. In section IV the proposed 
iteration scheme is employed for computing the correction to the zeroth-order approximation of 
the energy. The convergence of all integrals is demonstrated and the character of the singularity 
of the energy as a non-analytic function of the coupling constant is determined in the weak 
coupling limit. In addition, the details of all calculations are presented in appendices. 


II. MODEL DESCRIPTION 

Let us examine the Hamiltonian of the system consisting of a non-relativistic particle inter¬ 
acting with a scalar quantum field 

H = H 0 + H int , 

Hp = ^ a> k a l k a k , 

k 

Hi nt =-7= Y A k (e ikr a k + e~ ik r al) 

V2Q V v k ’ 

Here, we select the system of units in which m = 1 ,H = c= 1, the momentum operator p = 
-iV, normalization volume Q, vertex function A k = 1 / y/ZZyT, creation (annihilation) operators 
a' k (a k ) of the field mode with the frequency u> k = k = |fc|, and the coupling constant g. The real 
physical system, which is described via Hamiltonian (1) corresponds to an electron interacting 
with acoustic phonons in a continuous model of a crystal [17]. If we choose <x> k - 1 ,A k = l/k, 
and g = V87ra, operator (1) corresponds to the Frohlich Hamiltonian [23], which describes the 
interaction of an electron with optical phonons in a crystal, i.e., the so-called “polaron” problem 
[24-26], 

The total momentum operator 


(1) 

( 2 ) 

(3) 


P = ~ iV + E fca fe 3fc (4) 

k 

commutes with the Hamiltonian of the system (1) and consequently the eigenvalues E(P ) and 
eigenfunctions |'Tp) are defined as solutions of the following system of equations 

Hl'Pp) = E(PWp), (5) 

PITV) = P^p). (6) 


In the conventional perturbation expansion over the coupling constant in the zeroth-order ap¬ 
proximation the solution of the stationary Schrodinger equation with Hamiltonian (2) is simply 
determined and corresponds to the free particle with momentum p and Fock states of the phonon 
field with the set of occupation numbers {n kl , n kl ,...[ = {n k }: 


I'Ck,) = 2 a l a *u n *» = E nk]{,Jk]) * 

v i 2 Iff u 


E (Q HP,{n k })= L 


p - e kn k +E wfc " fc ’ j°=p+E feHfc - 


(7) 

( 8 ) 


Let us suppose that the system is in the ground state of the phonon field {«&} = 0, which leads 





4 


to the following eigenfunction and eigenvalue 


ivp(O) 

|T P,(> 


> = 


£ (0) (P, 0) = 


Vq 

p 2 

T’ 


|o> 

p — p. 


(9) 

( 10 ) 


The hrst non-vanishing correction to the system energy arises in the second order of perturba¬ 
tion theory (single-phonon intermediate transitions) and corresponds to the self-energy diagram, 
which defines the mass operator £(P) and is determined as 


Z(P) = A£ (2) (P,0) = 


1 


1 


a>k k 2 /2 - P • fc + (Ok 


— f 

167T 3 J 


dk 


16tt 3 J k[k 2 /2 -P k + k ]' 


( 11 ) 


In order to select bound state energy Eb = E a> (Q, 0) and effective mass m* of a particle we 
expand the energy in a series over P up to second order 


£ (0) (P, 0) + A E a \P, 0 )~E b + 


P 2 
2m* 


g 2 r dk 

16tt 3 J k[k 2 /2 + k] 



16 7T 3 


f 


dk 


k[k 2 /2 + k] 3 

( 12 ) 


(P ■ fc) 2 . 


The hrst integral in equation (12) logarithmically diverges, that is the bound state energy 
depends on the momentum cut-off K 


£ » = -^ ln 



(13) 


and becomes infinite when K —> oo, such that the correction to the energy is undefined in the 
framework of the perturbation theory for our model [27], 

At the same time, the corrected mass is well defined and equal to 


2-.! " 2 
m* 


g~ 
6n 2 ’ 


8 

6jt 2 


(14) 


In contrast to this in the polaron problem all integrals are convergent because they contain 
in the denominator the additional power of k. The corresponding quantities for the polaron 
problem read as [24—26] 


Eb - ~cr. 



(15) 


It is important to stress here that in our model, the interaction energy between particle and 
held is observable and consequently, the infinite energy (12) can not be included in the mass 
renormalization. Thus, we can conclude that the use of perturbation theory for two physically 
close quantum-held models leads to qualitatively different results. Therefore, a modihcation of 
the calculation method of subsequent corrections to the energy for our model is required and 
appears achievable. 


III. ITERATION SCHEME, BASIS CHOICE AND ZEROTH-ORDER 
APPROXIMATION OF THE SYSTEM’S ENERGY 


In order to build an iteration scheme not in the framework of perturbation theory we will 
employ the operator method (OM) for the solution of the Schrodinger equation, which was 
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introduced in reference [28] and its detailed explanation is given in the monograph [13, 14]. Let 
us quickly revise here the basics of this method. Suppose, the eigenvalues E M and eigenvectors 
IT^) with a set of quantum numbers /./ of the stationary Schrodinger equation need to be found: 

Hl'fy = fvny. (16) 

In contrast to perturbation theory, where the Hamiltonian H of the system is split into the 
zeroth-order approximation and perturbation parts, according to OM the total Hamiltonian is 
taken into account as is, while however, the state vector is probed via an approximate state: 

1%) * I 

which depends on a set of variational parameters a> II . Then, the exact solution can be represented 
as a series 


nv = + X (17) 

v+H 

Here we want to pay attention to the fact that for a given set of quantum numbers fj, the set u)^ 
is fixed. By plugging the expansion (17) into Schrodinger’s equation (16) and projecting on dif¬ 
ferent states |i/f /J (« /J )) and one obtains the equations for the energies E^ and coefficients 

r ■ 

'-'/iV 



1 + ^ C^ylfjy 

-1 





v^/j. 


& 

II 

F 

i 

i 

Hyn Ejjly^ + ^ ^ C^Hyy Eylyy) 


v+n+y 


tifJV — (t^/j(tUyj)|H|(/^y((U^)), /yy — (ffl I fv (^/i)) • 


( 18 ) 

(19) 


It is important to stress here that all matrix elements are calculated with the full Hamilto¬ 
nian of the system and the set of vectors | can be normalized, while not necessarily 

being mutually orthogonal. The system of equations (18), (19) is the exact representation of 
the Schrodinger equation. For the approximate solution of this system, in accordance with OM 
the following concept is used: the closer the zeroth-order approximation of the state vector is 
to the exact solution, the closer the matrix H )lv becomes to the diagonal one. Therefore, an 
iteration scheme for the solution of the system (18), (19) can be built, for which convergence 
is determined with the ratios of non-diagonal elements H^, to the diagonal ones H m in the 
representation of the state vectors |i/^(«/i)). A sufficiently detailed discussion of the conver¬ 
gence of the iteration scheme for different physical systems is given in the monograph [13], 
Consequently, we find the system of recurrent equations 


-t 


Ef = 

i+2^"% 

V±(l 


H + V' 

v&n 

> 

C % = 

E^-HyyY 

= 0; = 

L1 yn ^fj 1 yfi t '-'fjv V-^yv ^ 1 yvJ 

~-H m . 


( 20 ) 

( 21 ) 

( 22 ) 


As opposed to conventional perturbation theory, where the exact solution is defined as a sum 
of corrections of all orders, in OM the exact value of the energy of the system is given as a limit 
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of a sequence 


En = lim E^\ s = 0,1,.... 

S —>00 r 1 


( 23 ) 


In particular, for the first two iterations one can find 


P(l) _ p(0) _ TT . 


p(2) _ 


1 + 


Z 

v+n 


(.H, 


V[1 


E (0) I )I 

J-'H 1 VI2j 1 fiV 


E (0> - H 


+ 


Z 

V#/i 


(W, 




77 ( 0 ) 7 \U 

■ Ejjj 


E i0> - H 

1 - , H 1J v\ 


(24) 

(25) 


The last equation looks analogously to the second-order correction of perturbation theory, 
while the main difference is related to the denominators of equation (25), where the matrix 
element Hyy is calculated with the full Hamiltonian of the system, whereas the perturbation 
theory relations merely involve the diagonal element of the unperturbed Hamiltonian. As will 
be shown below, this is exactly the reason, which determines the convergence of integrals over 
intermediate states. 

Before proceeding with the application of the iteration scheme, let us discuss the choice of 
the parameters ,..., w'',...) in more detail. For this we note that the representation 

(18-19) is exactly equivalent to the Scrodinger equation (16), provided that the set of states 
{|<MK})» a complete one. It is evident that if the state vectors {|(A a ,({« a ,}))} coincide with the 
exact eigenstates | V F /J ) of the full Hamiltonian, the matrix // ;n , is a diagonal one, i.e. = E^d^ 
and the coefficients C f , v = 0. The eigenvalues E M are determined exactly and are independent of 
the set of parameters {n^}. Therefore, the relation 

d E„ 

—£=0, n — {1,2,...) (26) 

8m" 


holds identically. 

According to our initial assumption we choose the trial parameters { u >in the basis states 
|i/r A ,({« A ,})) such that they determine the best possible approximation for the exact solution IT^) 
in the chosen class of functions. This is equivalent to the supposition that the off-diagonal 
elements of the matrix H jn , are small numbers such that the ratios f/ ; ,„///,,,, are proportional to 
some effective small parameter e. Therefore, the zeroth-order approximation of the operator 
method is chosen as 

(H)) = C<®= 0, |^ #) ) = lft(H))). (27) 

However, the matrix H^ v contains small off-diagonal elements, which need to be taken into 
account. Hence, the subsequent approximations read 

OO 

E„ = e s £' s) ({^}), (28) 

S= 1 
OO 

C> f () = 2 e s C' s v >({^}), /if v. (29) 

s= 1 


As the left-hand side of equation (28) does not depend on the parameters {} it is natural to 
require that in each order in e the right-hand side also does not depend on (w ;( ): 


dE„ 
doj" 


= 0, 


« = { 0 , 1 ,...} 


(30) 
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for each w''. In the monograph [13] it was demonstrated that the recalculation of the parameters 
\ojfi ) in every order in e speeds up the convergence of the iteration scheme, however does not 
change the qualitative behaviour of the energy levels of the system. For this reason, in all 
calculations below we will fix the parameters {«,,} via the zeroth-order approximation: 


dE 


= 0. 


(31) 


In what follows we apply the above approach to the description of our model. In accordance 
with OM we choose a variational state vector, which incorporates the qualitative peculiarities 
of the system. From a physical point of view a held can be considered as a system of an infinite 
number of harmonic oscillators. Due to the interaction with a particle the equilibrium positions 
of these harmonic oscillators are modified. In a representation of creation and annihilation oper¬ 
ators the shift of equilibrium positions corresponds to the displacement of a classical component 
M;, on these operators [29], i.e. a' h —> a' k + u* k and a*, —> a*. + uk, such that we choose a basis of 
held oscillators consisting of coherent states. As a result a so-called localized state of a particle 
in the held of these classical components arises. This means that during its existence the particle 
becomes “dressed”, i.e. somewhat smeared out while still localized. Moreover, this “dressed” 
state should be an eigenstate of the total momentum operator P. since P commutes with H. 
Concluding, we formulate the following conditions for the state vectors: i) representation in the 
basis of coherent states; ii) imposing the localization of the particle state; iii) use of variational 
state vectors as eigenstates of P (4). 

In order to incorporate the hrst two conditions in the state vector, we choose it as the product 
of the square integrable wave function of a particle, localized near an arbitrary point R in space, 
and a coherent state of the held, analogous to the polaron problem [24, 30]: 


| t F(r, R )) = (f>{r - R)e \p 




,, J kR 

3i, - 


a fe) 


k 


| 0 >. 


(32) 


In the state (32) the classical component of the held Uk and the wave function 4>(r - R) can 
be considered as the variational parameters {o^} of OM. In accordance with the above described 
procedure of the choice of the parameters jw ; ,), the functional derivative over these parameters 
from the functional (Ttr, ff,)|H| v F(r. R)) should be equal to zero. This yields an equation for 
the classical components of the held iik and the wave function </>(r - R): 

[fflr, R)\mr, R))] = — 6 [0F(r, R))] = 0. (33) 

oiik o<p(r - R) 

By calculating the functional with the Hamiltonian (1) and corresponding derivatives one 
obtains the connection between the classical components of the held and the wave function of 
the particle: 

u k = - j 1 - f dr\cf>(r)\ 2 e- lkr . (34) 

l 

In the general case, the second equation in (33) leads to the integral equation for the function 
<pp(r). However, according to reference [13], the convergence of the iteration scheme of OM 
does not depend on the particular choice of variational parameters, under the condition that the 
approximate state vector takes into account qualitative characteristics of the system. Therefore, 
for the analytical investigation of the energy E ( k \P,g) we replace the exact numerical solution 
with a trial wave function, which depends on the single parameter A and is equal to 

A i 

<p(r) = —e ^ . 

7M 


( 35 ) 
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We notice that in the polaron problem the application of OM with the trial wave function 
(35) yields an accuracy of the order of 1% in the calculation of the bound state energy and the 
effective mass [30]. With this choice of wave function, we proceed to calculate the classical 
component of the field it k and the Fourier transform of the wave function (35), which will be 
required below: 


= - 


g 


1 


V2H VP 


f 


dr\4>(r)\~e 


2 -i k r 


g e « 2 

V2Q VP' 


tp k - ^ drtp(r)e lk r = 2 u 2 

J A 2 


2 A 2 . 


(36) 

(37) 


Furthermore, the states (32) are not the eigenstates of the total momentum operator P of the 
system, i.e. they are not translationary invariant. Moreover, these states are degenerate, as they 
do not depend on the localization point R of the particle in space. The choice of the correct 
linear combination of these states allows one to build a set of states which are not degenerate 
and are eigenstates of the total momentum operator P: 


l'P^ jnfe ) = — - 1 —j= J'dRcf> Pt (r - R) exp [i(P\ - kn k ) ■ R] exp j^(M g e lqR a' q - u q e lqR a q y^ \n k ), 

(38) 

P| v F l p > = P 1 |'f , p ) ). 

1 Pi,ilk' 11 Pi,rik' 

Here Q is the normalization volume and Pi the total momentum of the system, \n k ) are Fock 
field states with occupation number n k , (f>p l (r - R) is the wave function of the particle localized 
at point r — R and the classical component of the field u k is dehned via equation (34). The 
normalization constant for the state It'd 1 . ) is defined as 

Pl ■ 1 k 

\Np ak \ 2 = f dRidpt^pWnfp- fliy (Pi - fc, - Bl+j: ‘ |,, ‘ |2( ^ fc ' Rl - I, (2|K t [ 2 (cosfe ■ Ri- 1) + l). 

(39) 


In addition the set of states (38) forms a complete and an orthonormal basis in Hilbert space. 
The completeness of these states follows from the fact that they are eigenstates of a Hermitian 
operator P, with an explicit proof given in Appendix A. Concluding, the set of states (38) takes 
into account physical peculiarities of the system, forms the complete set of states in Hilbert 
space for arbitrary functions <pp(r ) and u k , which are an analog to the parameters {w ;j }, and, 
therefore, can be usable in the iteration scheme (20-22). 

The zeroth-order approximation for the ground state vector following above procedure then 
reads as 


IT^) = 


— f 

3 Vo J 


dR<pp(r - P)exp 


iP R+Z( 


u u e 


—i k-R 


a k - u k e 


,i k-R 


N P Vo 

whereas equations (24), (25) look like 

Ef + I,P M] C% {nk] {^ IHITV,,^,) 


7(2) _ 




(t-)ll 


^Pl,{n k ) 


Hp u {n k )-,Pi,{n k ) - E\ 


(0) 


HpUnkUPMik) = <^P„i» fe )|H|T'p 2>lfc) ), 4°' = ( v T , p|H|T' l p ). 


t/Wii 


a k ) 


U (L)\ 


10), (40) 


(41) 


( 42 ) 
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We want to emphasize once more that all matrix elements are calculated with the full Hamil¬ 
tonian of the system 


H = 


1 


\2\ 


pl - 2 Y a i a * k ■ p + Y a I a fc fc + Y wfca fe afc + E ( e ' fcrafc+e " fcr a 




/ / fe 


(43) 


Let us calculate the ground state energy of the system in this basis. The details of the 
calculations can be found in appendix C. The ground state energy reads accordingly 


Ef\P,g) = — -PQ + G + Ef(P) + E- m (P), 


(44) 


with 

Q = 2 fc|Mfe|2 f dRdr ^ r ^P (r - R)e^ R)+i{p - kyR - 

G = I ‘ V m ■ l\u m \ 2 \in\ 2 f drdRcp*(r)<p(r - i?yPR+*(RH 

Ef(P) = Y [ k + y) l Mfc l 2 f dRdr <p* p (r)ipp(r - R) e ' s,(Rmp - k)R - 
E m i(P) = ^p Y = f dRdr {tptr + R)4>p(r) + <p*p(r)(pp(r - R)) e *(*)*(P-*+*-->; 

o(ii) = 2lMfclV ifcfl -i); 

k 

\N p \ 2 = J dRdr <Pp(r)<p P (r - R)e" (R]+ ' PR . 

Actually, the iteration scheme (20), (21), (22) can be used for arbitrary coupling constants 
[13]. However, as was described above, in the framework of our model we are interested in the 
behavior of the ground state energy E ( ^ )] in the weak coupling limit. In this limit we can neglect 
the function 

d>(R) = Y\»™\ 2 {e- imR -l)~g 2 , 


in the exponent of all integrals in equation (44) as g <sc 1. 

First of all, we investigate the situation of a particle at rest, i.e. P — 0. In this case for 
the weak coupling limit the integrals in equation (44) can be expressed through the Fourier 
transforms of the wave function (f>{r ): 


dRidpf(p)<p(p- R l )e- ik Rl = J ~ dpf(p)e~ ikp JdRcf>(R)e ikR 


-fc = <Pl (45) 


J" dR\dp(p*(p)(j)(p - i?i) = |0 O | 2 = <p], (46) 
j " dRidpf(p)cp(p - Ri)e~ lk p = <p* k <p 0 = fikfio- 


( 47 ) 
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With the use of equations (44), (45-47) we can rewrite the energy of the ground state in a 
form 


4 0) (°’S) = 


l ym-l\u m \ 2 \ Ul \ 2 ^L 

2 ^i n 



2 g y Uk 4>k 

V2^ fc V*A>’ 


(48) 


which up to fourth order in g yields 


E?\0,g) = (A(- 4 + V2) V3 k + A 2 ) + 0(g 4 ). 


24n 2 

By minimizing the energy with respect to A one finds 

,(-4+ V2) 2 


E ( °\0,g) = -r 


32n 


/l= ^p( 4 - V2). 


(49) 


(50) 


In the weak-coupling limit it is also possible to obtain a renormalization for the mass of the 
particle. This is accomplished by expanding the energy (44) in a series over P up to second 
order. The details of the calculation can be found in appendix D. The result reads 


E?\P.g) 


7 ( 0 ) 


p2 


g 2 17 - V2 
9tt 2 21 


„( 0 ). 


1 + 


g 2 17 - V2 
9tt 2 21 


(51) 


From this equation we can conclude that the factor, which determines the corrected mass is 
half the one via the leading second-order term from perturbation theory, see e.g. equation (14). 


IV. SECOND ORDER ITERATION FOR THE ENERGY AND CONVERGENCE 

In the previous section we have found the energy of the ground state and the renormalized 
mass, which are proportional to the square of the coupling constant in zeroth-order approxima¬ 
tion. However, the correction to the energy coming from single-phonon intermediate transitions 
is of the same order with respect to the coupling constant. Consequently, its contribution should 
also be taken into account, thus requiring the calculation of the energy of the system in the 
second iteration (41). 

In order to calculate the second order iteration for the energy we notice (appendix E) that the 
matrix elements (T'p, ,{n k }\ 'Pp ) and ( l Tp l j„ fc ||H| l P ( p ) ), which are found in equations (41), (42) are 
proportional to the delta function of the total momentum of the system 6(P\ - P). Therefore, 
during the evaluation of the sum over P\ in equation (41) for the energy we have used the usual 
procedure [8]: one of the delta functions in its square was replaced through the normalization 
volume Q, and the integration over the remaining one yields P = P\, thus expressing the 
conservation of momentum. 

Firstly, we consider the case, when a particle is at rest, i. e. P = 0. The results, which are 
expressed through the Fourier components of the particle wave function in the weak coupling 
limit read: 


E (2 \0,g) = 

D 


(52) 
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where 


^ = 4 0) + X 


i 


r 1 w 

X 


2 


-Mfcrf’l, | p +k\- 


8 ( f>k<f>0 


u k<t> 2 k (y +£| + 


V 20 Vk 

g <Pk<P() 


- « k<P 2 k (g 2 lk + g 4 Jk - £/ 0) ) 


V2Q Vk 


+ Uk<P k (g~I k + g 4 Jkj - Efuktl 


'( 0 > 


— + ^| + ^ Ik + g Jk - E l 


(-D 


and 


B = 




(53) 


1 


(¥+*) 

_ g Mo 
Wn Vf 

- i<k(/> 2 k 

'g 2 I k +g 4 J k -E ^) 

«fc(<4 “ 0 q) 

0 

,1 +k , 

+ g 2 Ik + g 4 Jk ~ Ef ] 


(54) 


In equations (53- 54) we have introduced the following notations 


2 ™i«rj 2 4> 2 m+k = r44 n ; I( k = 

m 

z(f + = «> = 


k A 2 4 k - e 3 32 V6PlErf- 

^2 32^ k 


/F* 


T 2 v/foVe s 5 Erf ^ + 6^Erfi ^ 


(55) 

VP 




V2Q^ 


V ^(0™ +fc +< 

* ^ a / m 


v„- fc ) - V44 3) ; 4 3) = 


96tt 2 


v Erfi ii 


47t 


4 = fc-J^ + /® + /f; 

5 1 ' 2 a 2 


^ l ■ m\ui\ 2 \u m \ 2 <P l 


'l+m+k ~ g 4( P k Jk', Jk 


2 -l£ _ 1 

15 ,i 2 1 

e 5 ■ 


4(27t) 3 3 5 (1 + A|)3 


(56) 

(57) 

(58) 

(59) 


where Erf(x) = 2/ yfn ^ e~~dz and Erfi(x) = -iErf(ix) are the error function and the imaginary 
error functions, respectively. When we calculated the energy of the ground state, we dropped all 
terms with power in g larger than g 2 . Consequently, we can neglect the term g 4 J k in comparison 
with g 2 Ik , which can be confirmed by the direct numerical calculation of the integral. 

Prior to the numerical evaluation of the integrals (53) and (54), let us understand their struc¬ 
ture through the approximate analytical calculation. We investigate the behavior of the numer¬ 
ator and denominator of the quantities A and B. We start from breaking the integration region 
into two parts, namely [0, k () \ and [k () , 00 ). The value k {) will be fixed below. Let us work out the 
behavior of the quantity I k for small and large values of k. First of all we notice that g 2 / 0 gives 
exactly the ground state energy E (4)> . For small values of k, with the increase of k the value of 

g 2 I k - g 2 k 2 /(187T 2 ) + Ef\ i.e, it grows quadratically in absolute value, while being negative. 

Therefore, due to the presence of g 2 , this term is small in comparison with k 2 /2 + k for small 
values of k , so that, in the denominator of quantity A, the leading term is k 2 /2 + k. For large 

values of k, the quantity g I k exponentially grows as I k ~ e 3 * 2 /k and becomes the leading 
contribution in comparison with k 2 /2 + k, despite the higher power of g. 
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In analogy, we can analyze the numerator of the quantity A. For small values of A: we can 
neglect in every square bracket in equation (53) the large powers of g, i.e. terms with exponents 
larger than 1. Consequently, for small values of ft, the integrand within A looks like 


u k <p 2 k 

(f + ft 

t 2 

\ , g Ma 

) V2£2 V* 

0fc0o( 

f + ft 



(60) 


For large values of k, the numerator is exponentially decreasing, with the leading term being 
(—-^= ^^X-F^Mfc). This follows from the fact that ~ e and^fc ~ e Consequently, 
the integrand for large values of k can be presented as 


/ 8 0fc0O \/ Z7 1 

^ V2n V3E n ■ 

g 2 ( f>l Ik 


fV) 


( 61 ) 


Combining all together, we find that the quantity A can be approximately calculated as 

\2 


7(0) 


z 

k<ko 


\ Uk t (2 +k ) + V5n Vr) 

(f +*) 


■z 

k>k 0 


( 8 Z7' 

1 V2£2 yfk H 


(0) 


Ilk) 


<t>- k g 2l k 


(62) 


In this expression, both sums are well defined and remain finite. The sum over the region 
k > k[) is finite and convergent, while the ratio of numerator and denominator in the integrand is 

5ft 2 

exponentially decreasing as e . 

In equation (62) the point fto is determined as a solution of the equation 

ft 2 

-+k + g 2 I k -E ( £ ) = 0, (63) 

or employing the asymptotic behavior for the function I k (appendix F) 


ftg il 3 V67T e 3 -< 2 


and by finding the following logarithm 


In 


(f +ko)ko 


2 k~ Q 

= - 2 'lns' + 3 #’ 


(64) 


(65) 


with 


/l 3 V67T 

ci = - 7 —• 

48tt 2 

In the limit of extremely small g , we can build the solution of equation (65) via iterations, 
thus yielding 


ft 0 ~ T V3|lng|. 


( 66 ) 


The estimation of quantity B can be performed in a similar fashion and one finds 


B 


> + Z 


(' Uk % (2 +k ) + Uki< f > k 06 ) 


k<kt) 




z 

k>kn 


<p 2 k g 2Ik 


(67) 
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At first sight, it may appear that the quantity B features an infrared divergence, because the 
term /(0 2 (^ + &)) ~ 1/k 3 as k —> 0. However, this additional power of k in the 

denominator is cancelled through the difference <p 2 k - <p^ ~ k 2 / A 2 . The convergence at infinity is 

5k 2 

manifested with the exponential decrease of the integrand ~ e ’-j 2 . 

By plugging into equations (62) and (67) the values of 4>k and which are defined in 
equations (36- 37) and calculating the integrals (appendix F) we find the approximate analytical 
formula for the second iteration for the ground state energy 


(0) 


A* E 


B* 1 + 


8 1 * 


247T 2 


VfoErf 


y2^o 


\ \ 

3k 2 
__0 

A - Ae m 2 


g 2 'l 


2 V3tt 3 / 2 


-Erf 


V3k ( 

2A 


2n 2 


In I — + 11 + £' 


7 (0) 
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12 n 2 

i r x 

ml i 




ko\ 144V67T 


25An 


5 kf^ 


o 

12 A 2 


1 + 


e i 2 -i 2 , 


( 68 ) 

(69) 


tdt 

e 5 r . 


+ t/2 


Within the accuracy of the approximate formulas, we can set B ~ 1. Therefore, one finally 
obtains 


E (2 \0,g)*A. (70) 

The use of our simple analytical expressions allows to establish the behavior of the energy as a 
function of the coupling constant and consequently to determine the character of the singularity. 
In order to select the singularity, we investigate the limit 

lim £ (2) (0, g). 

In this limit, the value of ko logarithmically grows. Consequently, we can approximately set 
ko —> oo both in the expression in square brackets and in the last term of equation (68). This 
way, the square bracket becomes equal to the energy of the ground state (appendix F) and 
cancels E k \ The last term also does not contribute to the energy as being exponentially small. 
Consequently, the only term remains, which exactly determines the character of the singularity 
and is equal to 


<7i) 

ko ~ A y/3\ lng|. (72) 

We observe that this term exactly coincides with the result via perturbation theory, i.e. equation 
(13), however, here with a well specified “cut-off”. Moreover, the most contributions to the 
integral in the energy arise from the region k < ko and this is exactly the reason for the natural 
“cut-off”, which is determined self consistently and is directly related to the only parameter of 
the Hamiltonian, namely the coupling constant. Let us mention here that the corrections to the 
energy of the system (71) arise in the subsequent iteration and are related to the transitions into 
intermediate states with two phonons. These contributions are proportional to g 4 . 

In addition we note here that the absence of the ultraviolet divergence in the energy of the 
ground state in equations (52, 53, 54) is due to the fact that as in the zeroth-order approximation 
and in the second-order iteration in the resolvent [£^ 0) - H vv ]~ l of equation (25) the dressed 


12 V67T 5t 0 

- e ru 2 , 

5/br 
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Figure 1. (Color online) (a)) The dependence on the coupling constant of the ratio of the exact numerical 
evaluation to the approximate analytical formula of the second iteration for the energy. The value k 0 in 
the analytical approximation is equal to ko = A -y/3| lng|. (b)) The dependence of the imaginary part of the 
energy of the system on the coupling constant. The imaginary part corresponds to the finite lifetime of the 
state. 


wave functions (38) were used. This leads to the effective momentum cut-off k n (g), which is 
determined as the solution of equation (63). This cut-off is a function of the coupling constant 
and is not a phenomenological parameter, which needs to be introduced for the removal of the 
ultraviolet divergence. Moreover, as follows from equation (71), the energy of the ground state 
has a logarithmic singularity as g —> 0. It is clearly seen that this dependence can not be sorted 
out in the framework of perturbation theory, which yields a power series over the coupling 
constant g. 

In order to ensure that our interpretation is correct, we have evaluated the integrals numeri¬ 
cally and have found in the limit of extremely small g the ratio of the results via exact numerical 
versus analytical evaluations. This ratio is almost constant and is approximately equal to one, 
as presented in Figure 1. Therefore, we can conclude that the main reason why conventional 
perturbation theory fails is related to the fact that the energy of the system is a non-analytical 
function of the coupling constant and consequently can not be expanded in a series over g near 
a singular point. 

The second interesting consequence of the numerical evaluation of the integral is related to 
the fact that the energy of the system contains a small imaginary part, which means that the 
state has a finite lifetime and is quasi-stationary. To prove this, we have calculated the transition 
probability to the state | v Fp 1- i fc ) for the case when a particle is at rest, i.e. 




O k 2 

u k ( t’ k ( 2 + ^ 

) + ^n‘ t $ L+ u ^lS 2 Ik - Efu k (t> 2 

2n\k+l+g 2 r k \ 


^i 


The result of evaluation is presented in Figure 1. As can be seen from the figure the two 
curves coincide exactly. This can be interpreted via the diagram technique [8], The second order 
iteration for the energy of the particle can be presented via the diagram depicted in Figure 2. 
If the diagram is split by the dashed line, the imaginary part will correspond to the transition 
probability to the state | v Fp 1 j fc ). 

Here, we need to stress that contrary to standard perturbation theory in our formulation the 
conservation of energy is governed not by the free Hamiltonian Ho, but through the expectation 
value of the total Hamiltonian F/p i fc: p, i fc . Therefore, for certain values of k and for certain 
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Figure 2. (Color online) Feynman diagram of the process. If the diagram is split by the dashed line, 
the imaginary part will correspond to the transition probability to the state |'Fp li i fc ). However, in the 
intermediate states the resolvent of the full Hamiltonian instead of a free one is used. 


coupling constants g this energy level might appear to be below the energy E ( ®\ featuring a 
so called quasi-intersection of energy levels. If the transition probability to the state | v Fp,,i fe ) 
were large, the description of the system with the state vectors (38) would not be applicable 
and the reconstruction of the states would need to be performed, which takes into account the 
degeneracy between the energies E^ ] and lip, i fc; p, i fc . In our case, however, the transition prob¬ 
ability is small and consequently the description with a complex energy, with small imaginary 
part, is valid, in analogy to the theory of a natural line width of the atomic states or anharmonic 
oscillator p 2 /2 + x 2 /2 - px 4 , with p > 0. 

In order to conclude our formulation we have calculated the renormalized mass in the second 
iteration. In terms of the introduced abbreviations the second order iteration for the particle 
energy can be written as 


E {2 \P,g) = 


2 


+ E ( °\P,g)+A P 


Bp 


(74) 


with 


5(0) 


(p, 8 ) = -p-Y a 


mkirr.r-^ + 
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,2'l'P-m 


The quantity Ap reads as 
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k 
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^ p2 Ik 2 \ 

Ap = ~ u k (fip-k ~ ftp) + + k — P ■ kj Ukfip-k + 
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0p 

' V2f2 
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V2f2 
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(75) 


+ 8 u k 4>~ P - k Ip-k - E k (P,g)u k (pp 
~b 8~ u k<t>p~ k Ip-k 


fp-Jp 


n — 1 


-+k-P-k\+g l I P - k -Ef\P,g) 


(76) 


and the quantity Bp is equal to 


Bp -1 + z 


E L \P,g)u k <pp_ k [I 0 +k P k\u k (j)p_ k + 


8 0P-fc0P 
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To proceed, we again break the limit of the integration into two parts, however, now we 
know that the main contribution to the energy of the system comes from the region [0, Atq]. In 
this region, we again drop all terms, with power of g larger than one. We recall here that the 
classical component of the field u & is proportional to g and the energy E ( ^(P,g) ~ g 2 . In 
addition, the limit P <K 1 is considered. Moreover, as the quantity Bp, after expansion over 
momentum P, will have a form Bp = 1 -g 2 F(P 2 ), we can set P - 0 in B p, in order to preserve 
the same accuracy. 

In this approximation, the quantity Ap takes the form 


Ap - 


P 2 vn f (*P-k~*p) 


z 


Uk- 


k<k 0 


bp-kfip 


4>p-k , 8 l k 

Ilk -1-;= —— ( — + k - P ■ k 


V 20 Vfe \ 2 
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k<kn 


4>P-k 


8 1 


- 1 - k - P ■ k Mfc ,_ 

2 I <Pp V20 Vk 
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— + k - P ■ k 
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(78) 


and 


B P = B 0 = 1 - 


( u ^{^ + k ) + Wc 


k<k 0 


8 <ph 

V 5 n 00 Vk) 
<Pl(T+k) 


(79) 


If the definition of E { ^\P,g) (75), together with equations (78) and (79), is used in equation 
(74), the second iteration for the energy of a moving particle is obtained 

E {2 \P,g) = E (2 \0,g) - — y ' fc) 

S 20 ZZ k(k 2 /2 + k) 2 
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(80) 


or after simplification, taking into account the fact that the sum in the denominator is propor¬ 
tional to g 2 , one finally obtains 


E (2 \P,g) = E^(0,g) + 
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11 

20 


Z 

k<k 0 


(.P ■ k ) 2 
k(k 2 /2 + k) 2 ' 


From here, we see that the second iteration for the renormalized mass 

„2 
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( 2 )* 
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8~ 

6ir 2 


(81) 


(82) 


coincides with the one via perturbation theory. 


V. CONCLUSION 

In current methods of renormalization in QFT, the momentum cut-off plays an important role 
[31], which in fact is an additional and undefined parameter of the theory. Usually, the inclusion 
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of such parameter for a concrete model is justified with the argument that the theory becomes 
incorrect on a small scale, where a more general theory must be used instead. For example, 
in the case of QED it is widely accepted that on a small scale the Standard Model, with its 
own characteristic length, should be rather used. However, in the Standard Model, as in its 
possible generalizations, for the renormalization of perturbation theory the cut-off is required. 
Consequently, we come to the requirement of the inclusion of some “fundamental length” or 
unobservable parameter of any QFT. 

However, the Frohlich Hamiltonian demonstrates the absence of a cut-off in the polaron the¬ 
ory. In this QFT all corrections are determined through convergent integrals and, consequently, 
the cut-off is not required. Here we considered a more general QFT than the one associated 
with the polaron problem, for which standard perturbation theory gives rise to divergences. The 
main result of the present work consists in the construction of a calculation scheme for this 
more general QFT that only leads to convergent integrals. In addition to that, the regulariza¬ 
tion of all integrals is related to the effective-cut-off momentum, which is defined through the 
parameters of the system itself. Moreover, the divergences of standard perturbation theory are 
explained through the energy being a non-analytical function of the coupling constant, of a form 
ln( yj | lng|/2 + 1), around zero, and, therefore, can not be represented as a power series around 
this singular point. It is also important that the character of the singularity, defined in equation 
(71) in the weak coupling limit does not depend on the particular choice of the wave functions 
4>p(r) of the zeroth-oder approximation. 

From a formal point of view, the convergence of all integrals is explained as follows: i) the use 
of the decomposition (17), i.e. the special state vectors, which are the product of the wave func¬ 
tion of a localized particle and a coherent state of the held and ii) the calculation of the energy of 
the system with the iteration scheme (25), in which the resolvent of the operator [H^ - £,’ (0| ] _l 
contains the matrix elements of the full Hamiltonian of the system. In standard perturbation 
theory the Hamiltonian of non-interacting fields is used in the analogous expressions. 

From a physical point of view the argument i) corresponds to avoiding an adiabatic switch off 
of the interaction. This means that a particle during its existence time is considered as “dressed”, 
i.e. to be in a localized state which is created due to the interaction between the particle and the 
field. The argument ii) leads to the “cutting” of all integrals for a large momentum due to the 
reconstruction of a localized state in intermediate states, caused by the quasi-intersection of the 
ground and the single-phonon states. 

Our approach should not be considered and does not pretend to be the full solution of the 
renormalization problem in QFT, specifically, because of the use of a simple, non-relativistically 
covariant model. Nevertheless, it demonstrates an alternative, succeeding without introducing 
any phenomenological momentum cut-off. 
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APPENDIX A: PROOF THAT THE STATES (38) ARE EIGENSTATES OF THE 
TOTAL MOMENTUM OPERATOR 

In this appendix we present an explicit proof that the total momentum operator 
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commutes with the Hamiltonian 


H = -ii + 2 («"' "*« + e ~““ r <) 


fin „ 

of the system and that the states (38) are eigenstates of P 


) = Pi |T® >, 

1 Pl.Hfc' 11 Pl./ifc'’ 


consequently forming a complete set in the Hilbert space. 
Let us begin with the commutator: 


(83) 


[H,P] = 
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V 2 D k 




a k +e- ikr a{),-iV 




Y J A k {e' kr a k+ e-' k -a[),Y J ^W 


V2Q^ 


Y J A k (a k e ikr (-k) + a {e~ ikr k) 


+ ^=Y J A k (a k e' k r k + ale-' k r (-kj) = 0. 


(84) 


which was to be proven. 

Now we will demonstrate that the relation (83) holds. Also, we will introduce the notations 


i = exp |^(Mq<? ' qR a l q - u q e iqR a q )j , with 

(85) 

D\R)D(R) = D(R)D ] (R ) = 1, 

(86) 
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(87) 

D\R)alD(R) = a; + u* k e ikR , 
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= DCR) Z ^ e ~ iqR < + »y qR a q ) 

(89) 


Consequently, with the help of equations (85-89) we may write 
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By noticing that -iV r (^p,(r - R)) = iVp(^p,(r - R)) and transforming = 
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D(R)D\R) qa' q a q D(R) one obtains 
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P|'F® ) ) = 

1 Pi ,«fc' 


N Pl,n k VQ 
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N Pl ,n k VQ 


— f 

„ Vo j 


N Pl ,„ k Vo 
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x X g ( a q + ll q e ' q R ) ( a <? + u q e~ iq R ) | n k ). 


(91) 


The first term in equation (91) vanishes due to the square-integrability of the function <p Pl (r - 
R). The derivative in the second term is equal to 


(iV fi ) [exp {i(Pi - kn k ) • R) D(P)] = -(Pi - kn k ) exp {i(Pi - kn k ) ■ R) D(P) 

+ exp {i(P[ - kn k ) ■ R] D(R) ^ q (u q e~' qR a\ + u* q e lq R a q ) 

(92) 


and, therefore, the terms which are not proportional to P\ in equation (92) cancel the last term 
in equation (91). As a result, equation (91) transforms into 


■n.Hfc' 


(93) 


which was to be proven. 

According to reference [32], the eigenstates of a Hermitian operator form a complete and 
orthogonal set of functions in the Hilbert space. As the functions (38) are eigenstates of the 
Hermitian operator P, they form a complete orthogonal set for arbitrary generalized parameters 
d)p ] (r - R) and u k . 


APPENDIX B: MATRIX ELEMENTS CALCULATION 


In all subsequent calculations, the matrix elements of a type 
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need to be evaluated. By using the identities 
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equation (94) can be transformed into the form 

exp(E m I u m \ 2 (e-' m {R - R,) - 1)) 

x <0|(aj - Uje~ lj R + Uje~' j R ) ni ^ /(at + uie~' l R , a]' + u\e ll R ) 

l 

x (a[ + u* k e' k R - M ;e ,fc R )' ifc |0). (97) 

The evaluation of equation (97) is performed in the usual manner, i.e, by noticing that a|0) = 
(Ola 1 = 0 and the vacuum average is not equal to zero only if the number of creation operators 
is equal to the one of annihilation operators and is an even number. 


APPENDIX C: GROUND STATE ENERGY 


According to equation (42) of the manuscript, the ground state energy is defined as 
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The normalization constant Np is found from the condition 

i'Vp’l'Vp) = 1 . ( 101 ) 

In order to evaluate equation (101), we use equation (97), in which n k = rij = 0 and /(a*, aj) = 
6ix). This gives immediately the result 
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( 102 ) 


The expectation value of the energy is performed in exactly the same way. First of all, the 
matrix elements of the field states are calculated with the help of equation (97). For example, if 
the function / is chosen as / = Za^a;, n k = tij = 0, we immediately find 
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Then by carrying out the change of variables R - R' — R\ and r - R' = p, we obtain 

Q = ITT"i2 X fc l Mfc l 2 f dRdr ^( r - ^)e° ( ' R)+i( ' P - fc) jR , (104) 
ItVpr ^ J 

0(R) = ^\u k \ 2 (e- ikR -l). 


All other matrix elements are evaluated in exactly the same fashion. Consequently, we obtain 
expression (44) of the manuscript: 
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We notice here one more time that the Fourier component of the function reads 


P- -iV 

<t>(r) = —e - . 
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(106) 


and the classical component of the field look like 
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In order to calculate the energy, we firstly neglect the function 
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in equation (105). The remaining quantities can be rewritten employing the Fourier transform 
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of the function (Hr). As for example 
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( 112 ) 


and by plugging equations (110-112) into equation (105), one finds 


(0,g) = ~ J] rn ■ l\u m \ 2 \ Ul \ 2 

1 ml *0 


(113) 
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By insertion of the definitions of the classical component of the field u k and the Fourier 
transform of the function (f> k , defined in equations (107) and (108), we find 
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(114) 


where a = 0.736559. 

To find d, we minimize the energy, which results in the equation 
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from here we find 
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and by plugging d in equation (114) 
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In the weak coupling limit we find the renormalized mass in zeroth-order approximation. For 
this purpose, we rewrite the energy through Fourier components for the case P + 0. This yields 


p 2 

Ef\P,g) = — -P Q + G + Ef(P) + P int (P); 018) 

Q ~ ^Yj.1 7 ^ (119) 

G = \]TTP (120) 

2 1^1 w 

Ef(P) = Wf 2 (* + I) (121 > 

Eintl-P) = . g |2 V (0P0P-fc + fipfip+k)', (122) 

|!VpI 2 V V2kQ 

|Ap| 2 = 0-p- (123) 


From here we can immediately conclude that the quantity G yields only a correction of the 
order of g 4 and can be neglected. Let us expand the Fourier transform of the function (Mr) into 
Taylor series over P up to second-order 
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(124) 

(125) 

(126) 


By plugging equations (124-126) into equations (119-123) and taking into account only the 
terms of the order of P 2 , we find for the vector Q 
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and for the field energy 
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In this expression, in round brackets the first and the last terms cancel each other after decom- 
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position of the normalization constant in the Taylor series in P. The result reads as 


Ef = S f (0) + -2 ^ I fc + — I \u k \ 2 e i 2 (P ■ k) 
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The remaining energy is calculated in exactly the same way 
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In a full analogy to the field energy Ef, the first and the last terms are cancelled. The remaining 
terms are 
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(131) 


By combining all results together, we find the equation for the total energy of the system with 
a renormalized mass, in the zeroth-order approximation: 


Sf (P,g) = Ef\0,g) + — 
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or by plugging in for A according to equation (116) we finally obtain 
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Concluding, the renormalized mass is equal to 


*(0) 1 . 
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g - 
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(132) 
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(134) 


APPENDIX E: CALCULATION OF MATRIX ELEMENTS IN THE SECOND 
ITERATION FOR THE ENERGY OF THE SYSTEM 

The calculation of the second iteration of the energy of the system requires the evaluation 
of the transition matrix elements < v P ( p ) |H| l Pp 1 , { „ fc} ), <'T < p ) |'Tp 1 ,i„ te }> and <Yp li{ „ fc) |H|T'p 1 , { „ fc j> from 
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the full Hamiltonian of the system, equation (100), with the function 
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The normalization constant in equation (135) is calculated with the help of equation (97) and 
has the form 

llVp.iJ 2 = J di^ld^(p)^(p-i^l)e ^fl - h > fll+ 2•. | “‘ |2 ^ ,, • ,,, - 1 ) 
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The calculation of the transition matrix element is performed with the help of equation (97): 
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We further obtain the cover integral 
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and the expectation value of the Hamiltonian 
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(139) 


Equations (137-139) are valid for arbitrary coupling constants. However, in the weak cou¬ 
pling limit they are significantly simplified as they can be expressed via Fourier transforms. 
Another significant simplification comes from the fact that the action of one field mode on 
the system is inversely proportional to the square root of the normalization volume O, that is 
Mfc ~ 1 / Vo. Consequently, such terms can be kept only within the sum. Within this approxi¬ 
mation, equations (137-139) can be rewritten as 


/in , vp( L)\_(2^(P-Pi)« fc (0p_ fc 


and 




O 


<l K _ (2t:) 3 6(P -Pi) 1 

O < Ppfip-k 

P 2 

1 .. / ,2 .2 


(140) 


»p<pp-k 


^ UkifiP-k ^p) + ( ^ ^ P ■ k ) lt k Cpp_ k + 


g <t>P-k<t>P 

V20 Vk 


^ P-m-k 


- u k (Pi - k) ^ 

m 

+ Uk 2 (y + + V Zi 1 ' 

m ' ' Lm 


g VH M 

+ 


V 2 Q tV V^ 


Z ,— <P P-k((p P-k-m + (pP-k+m) M k <PpE 

\ m 


(141) 


as well as 
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In equations (140-142) we have used expressions for the normalization constants 


Np - (ftp, Np li \ k = 


( 143 ) 


and sorted out the energy of the zeroth-order approximation 
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APPENDIX F: SECOND ITERATION FOR THE ENERGY OF THE SYSTEM WITH 

PARTICLE AT REST 


For the evaluation of the second iteration for the particle energy the knowledge of the behavior 
of different terms in expressions (140-142) is required. In order to determine those, we will carry 
out the summations over m. For the first sum we can write 
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and for the third 
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In the above expressions Erf(x) = 2/ \fn f 0 ' e z dz and Erfi(x) = -iErf(ix) are the error func¬ 
tion and the imaginary error functions, respectively. 

Consequently, we can introduce the abbreviations, which where used in equations (55-59) of 
the manuscript, namely 
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After the determination of the asymptotic behavior of the different terms, we can find the 
second iteration for the energy of the system 
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For the numerator we have 
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Further explicit calculations yield 
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